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GAS FLOW WITH STRAIGHT TRANSITION LINE* 
By L. V. Ovsianhikov 


On the basis of the solutions obtained by S. A. Chaplygin 
(reference l) , an investigation was made of the limiting case of 
a gas flow when the constant pressure in the surrounding medium 
is exactly equal to the. critical pressure for the given initial 
state of the gas> the results are presented herein. For a jet 
flowing out of the opening in a vessel with plane walls, it is 
shovm that equalization of the flow in the jet is attained at a 
finite distance from the start of the free jet, the line of transi- 
tion being a strai^t line . 


1. According to Chaplygin (reference l), every problem on the 
determination of the subsonic flow that is satisfied by some condi- 
tions reduces to the solution of the system 


^ _ 2T ^ 

Se - bT 

g-T 

^ ■ 2aT(l-T)W 55 

or the equivalent system 


> ' 


^ ^ 2T ^ 

' (1-'t)P 

^ ^ _ 2aT(l-T)P-H ^ 
B<iP g-T S'|r 




> 

J 


( 1 . 1 ) 


( 1 . 2 ) 


for the corresponding "boundary conditions. 

In these equations, = ^(6 ,t) and = 9(0 , t) are the stream 
function and the velocity potentia!., respectively^ t = where 


*"0b Odhom Gazovom Techenii s Priamoi Liniei Perekhoda." 
Prikladnaya Matematika i Mekhanika. Vol. XIII, 1949, pp. 537-542. 




2 


NACA TM 1295 


V is the magnitude of the velocity vector and is the maximum 

value of V that satisfies the given initial state of the gas; and 
9 is the angle hetveen the velocity vector and the x-axis. Further- 
more, p = i/(k- 1) is the adiabatic power of the fxmction of the 
density on the ten^ieratxire, where for brevity, 

a = ^ 

23+1 K+1 

so that the value T = a corresponds to the critical velocity 



In equations (1.3) and (1.4), Q is the relative quantity of 
flow in the gas jet and C is an arbitrary constant depending on 
the choice of the origin from which the values of. are computed. 
The magnitudes Zjj, Zj^q, and Xj^ are defined by the equations 
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Zn = 2n('^) = 
ZnO = ^n('^o) 


^ = 


Xn('T) = 1 


^ T yn’(*^) _ T Zn'('^) 




(1.5) 


(n — 1^2^ • • •) y 

y^ere yni"^) is that solution of the hypergeometric equation 
t(1-t) Yn" + 1 (5-2n-l)T Jy^' + pn (2n+l) y^ = 0 

which is obtained for T = Qj Tq denotes the value of ,t cor- 
responding to Po* ■ 


In Chaplyginis investigation, it was assumed that t 
this ease, equations (1.3) and (1.4) converge everywhere in the region 
of flow and give the propqsed solution of the problem. The case 
where Tq = a is cpnsiciered in the fpllowing development; 

It is recalled that the functions Zn(T) and xn(i) for 
satisfy the following inequality pf Chaplygin: 



zn :^o 


(1.6) 


/ d-T . qT / a-T 

y Mi-r) „i/5 a-,} ” V 9-a-n 

5 ' d-'') 

q =^\/232 (i+ 2 p) = constant 

2. It shall be shown that the potential pf the velocity cp given 
by equation (1.4) fpr T -hi^Tq 


(a) increases ^thput liirdt if 

(b) remains finite if T o = 
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It shall he assumed that the limit T — >Tq is effected by 
moving along some streamline. The following notation is introduced: 


4 ^ t(1-t)^P 


OO 

P(0>'T>To) = V cos 2n0 

n Zj^Q 

n=l 


It is readily seen that for 0 


Tq ^ a, then 0 -c ^ *=C1 . 


For the potential at the points of the x-axis ’ (0=0) , the 
following expression is obtained from equation (1.4): 




1 r drr 

t(i-t 


)P (I-t)PL 


- 1 + P.(0,T;To) 


( 2 . 1 ) 


On the. basis of the inequalities of equations (1.6). and (1.7), 
for the magnitude of P(0 ,TjTq), 


n=l n=l ... n=l 


n=l 

( 2 . 2 ) 


zn - 

^nO 


Assertion (a) now follows in an obvious manner from the second 
of inequalities (2.2) because zJz^q—^1 for t — »Tq and the sum 

of the series with the general term z^/nz^Q therefore increases 
without limit; whereas the coefficient preceding the term, because 
Tq ^ a, approaches a positive limit. • 


From assertion (b), it is seen that 


E 

n=l 


i = _ log(l-^) = - log 

n . • ■ 


Tg-T)^p 

7.o(1-To)2P_ 



NACA TI.I -1295 


5 


"by virtue of which, for Tq = a as T— »a in equation (2.2) on the 
left, the first component approaches zero and the second component 
approaches a finite magnitude equal to 


OO 



The same. reasoning, together with the first of inequalities 
(2.2) is repeated for P(0,Tja) for 0^0, so that the assertion 
is completely proven. 

3. From the previously proven houndness of the velocity potential 
it follows that at a finite distance from the opening BB', the jet 
is intersected hy a certain line L along -which T = a, so that the 
velocity is equal to the velocity of sound. It will now he sho-wn that 
at all points of L, 0 = 0. 

First, it will he observed that along any fixed streamline, 

0 varies roonotonically, as is true of the boundary streamline, and 
the transformation (0 ,t) — »>(<?, ^ i^) is a single sheet transformation 
at every interior point of the flow region. Next, hy fixing some 

value T|f = i|r for which i Q, the limits is approached as 

and 0 q = lim ,0 is set for T— »a and i|r = i|r; this limit 
. exists because of the finiteness of 0 ana the monoticity of its 
variation. 


Approaching this limit in equation (1.3) for i|f = yields 



that is, in all cases for 0^ 4 0, |'|| = ^ Q, -which contradicts the 

t I 1 ^ ^ 

assumption that Hence, the equality 0=0 must hold. 

It VTill now be sho-vm that the line L is straight: Along L, 

9 = constantj this result, is readily obtained .if the previously con- 
sidered transition to the limit is carried out in equation (1.4). 
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It is sufficient to note that every displacement in the plane xy 
is connected with the corresponding displacement in the plane cpijf hy 
the relation 


dx = ^ dq. - dt 

V ^ V(1-T) 


(3.i) 


from which it follows that in a displacement along L in the 
xy -plane, dx = 0, because in a displacement along L in the 
cp\|f -plane, according to what was previously proven, 6=0 and 
dq) = 0. The line L in the xy -plane is thus a straight line 
perpendicular to the x-axis. 


The equalization of the jet occurs along the line L. Behind 
this line the jet becomes uniform, flowing with constant velocity 
everywhere equal to the velocity of soixnd. 

The distance of the line L from the edge of the opening will 
be counted. Along the boundary of the jet, ijr = constant and 
T = a = constant, so that at the points of this boundary, equa- 
tion (3»l) assumes the form 

cos 9 Sq) 
dx = d6 

cr 

Substituting the expression for qp (equation (1.4)) taken 
for T = a yields 




OP 


(a) sin 2n6 cos 0 d0 


(3.2) 


Inasmuch as 


2 sin 2n0 cps 0 = sin (2n+l)0 + sin (2n-l)0 

integrating equation (3.2) from x = xg, 0='|’rt, to x = x^, 6 
yields 


= 


rtV„„(l-q)P 4n^^l 

V-*- 1 


Xn(o-) 



NACA TM 1295 


7 


If h denotes the width of the jet where it is uniform, then 
Q = V(sj.(l-a)P h, so that the required distance is finally obtained 


in the form 


h 


1 

It 


OB 



4n 

4n2-l 


Xn(o-) 


(3.3) 


4 . If any two streamlines of the obtained flow are teiken as the 
walls of a certain nozzle, the subsonic flow within this nozzle will 
be determined by equations • (1.3) and (1.4). This flow becom.es uniform 
at a certain distance and has a straight transition line. This fact 
corresponds entirely to the result obtained by S. A. Christ ianovich 
in his investigation (reference 2), where a general device is given 
for the construction of a Laval nozzle with straight transition line. 


The first derivatives of 6,T vxith respect to 9,\lr are 
estimated near the line L inasmuch as the behavior determines the 
possibility, of continuing the obtained flow across L as a super- 
sonic flow. The investigation of this complicated problem has been 
started and, it is hoped, will be presented in a forthcoming report. 
Remarks herein will be restricted to the following: 


In the first place, from the previously noted properties of the 
line L, it follows that the derivatives 5 i/^\|f, and B0/Sqp 

are equal to zero on L. 


In the second place, the derivative dx/Scp is evaluated at the 
points of the x-axis. It is sufficient to consider the derivative 

inasmuch as, on the x-axis, because 0 = 0 and i|r = constant. 


5cp Bcp/ST 

Substituting equation (1.3) for ^ in the second of equa- 
tions (l.l) for T = a and setting Zj^^^ = Zj^(a) yield 




Q(g-i) 

2itaT(l-T)0+l 


•E 


^n 

cos 2n0 

2na 
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or, on the x-axis where 9=0 

^ = Q(g-T) 

jtaT(l-T)P+l 

For estimating the value of the expression in brackets, which 
is denoted by S, the last of relations (1.5) and the inequality of 
Chaplygin (equation (1.7)) yield 

• 



n,/IEII + _l_n2/3p,Alog Zn(T):^“^/I^ 
^ Vo'(l-'f) 1-7 17 7 ^ 0,(1- 


(4.3) 


Integrating equation (4.3) from to T = a yields 



pa 

pa 

1 / ^-7 dT + an2/3 

_^>log n 


Va(l-7) 

1-7 Zn 

V‘ 

>T . J 

7 J 

T ’ 


dT 


Raising the upper and lowering the lower limits and carrying 
out the integration give 


2n 


3T]^ ^a(l-a) 


(a-T)3/2 + qnV^(g-T) ^ (a-T)V2 


1-a 


‘-n oa' 




If a-T = z(0<^z^a-T]^), where z— as T — *a and 

2 


(4.4) 


r = 


3Ti(^a(l-a) 


>*0 




5 = 


->►0 


3aV^ 

equation (4.4) assumes the form 


r nz^f^ + q-n^/^zs^log flE;^6nz^/^ 

X rr 
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Then the inequalities 


exp( -rnz V^-q, n2/^z ),^_^^exp( -SnzV 2) 

rr 


(4.5) 


“net 


are obtained . 

Because of equation (4.5), the following inequeilities are 
obtained for 3: 

OC> , QD 


exp(-Y'nz'^/2-q^n2/3z)5;^S^^ + exp(-6nzV2) = S2 
n=l n=l 


The sum S 2 is readily cornputedj 


2 2 


exp(6z^/^)-l ^ 


(4.6) 

(4.7) 


In order to estimate the svxa S]_, ' 


Un = 

exp(-qQ^n2/3z) 


^n = 

exp( -vnz3/2) 


Sn* = 


\ 



k=0 



= 



exp f-q^n^/Sz-ynz^/^ 


n=0 

n=C 

W V 

Applying the 

transformation 

of Abel to yields 


m-1 



' 

% =YZ 

(%"%+l) ^n’ %i^' 


n=0 




( 4 . 8 ) 
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For 03 ^z^(id+ 1) o), where oi = r"^/^(log 2)“/^, 

for the Slim Sn' 

n 


Sn’ = ^ exp(-rkz3/.3)2>H^ 
k=0 


(4.9) 


The differences Un-Un+i are evaluated hy the formular of finite 
increments . Therefore 

' ^ “ "n+1 = e:q) [^-qj^n2/3^ _ exp q3^(n+l)2/3 ^ 

— (n+^)"^/^ exp zj 

3 


Replacing ® hy 1 yields the inequality 

(n+l)-l/3 exp[-q2^(n+l)2/3 ^ 

Comhining equations (4.9), (4.10), and (4.11) yields 
m-1 • m 


(4.10) 


itA -L. 

gm j qjz(n+l)2/3 exp ^qT^(n+l)2/3 zj = | ^ ' q2,n2/3 z exp jj-q^^nS/S ^ 

n=0 n=l 


n=0 
The function 


-(x,y) = xy2/3 exp(-xy2/3) 


(4.11) 


is considered for x>*0* The derivative \j±th respect to y 

varies as ^ 

-3/2 


dy 


> 0 

for 

O^y-c X 

= 0 

for 

y = X 

0 

for 

y> X 


For fixed • x>0, the function f(x,y) increases at first from 

0 to e-1, then decreases and, for y — »oo, approaches zero. Hence, 
for any x>»c, the inequality 
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1 + E n2/5 


n=l 


is obtained. 


X exp(-n 2/3x); 


m 

X5I-2/3 exp(-xy2/3) dy - 


(4.12) 


The value of the integral g is obtained by applying the 
substitution xj'’^/^ = t. ^en 

'X 


g = 


2^2 


-^3/2 g-t (X=m2/2x) , ; . (4.13) 


Substituting equation (4.13) in eauation (4.12), setting 
X = q-j^z, and noting that yield 


q-.n2/3 z exg{-q;\n^/^z)z 


2qi3/2 


n=l 

hence 


q.103 


0 


t2/3 Q-t d.t 


-3/2 _ 


' Si = - 1 -J- I + 


2qi3/2 


q-,0) 


t2/3 g-t 


dt 


0 


z-3/2 (4.14) 


On the basis of inequalities (4.6), (4.7), and (4.14), it may 
be concluded that there are two positive constants 6i and such 

that the inequalities 


6 iz 3/2 


^ <S^ \ 

y-, 


/2 


(4.15) 


hold over the entire interval O^z^a - T^. Returning to the 
variable t and comparing equation (’4.2) with equation (4.15) gives 
the following result: Two positive constants 63 and r2(°2^^"2^ 

exist such that for emy t in the interval Ti<;:-T«cc,, 


1 op 


62 TaV ^ 


(4.16) 
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From equations (4.1) emd (4.18) there follows finally 

62 ya-T>^>T 2 /y/(a-T’) (4.17) 

which is proven for the points of the x-axis. 

* / 

In a similar manner, it may be shown that on the x-axis the 
second derivative remains finite as T — »a. 

. Translated by S. Reiss 
National Advisory Coinmittee 
for Aeronautics 
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